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Abstract 

The existence of a dynamo effect in a simplified magnetohydrodynamic model of turbulence is con- 
sidered when the magnetic Prandtl number approaches zero or infinity. The magnetic field is interacting 
with an incompressible Kraichnan-Kasantzev model velocity field augmented with a viscous scale cut- 
off. An approximate system of equations in the different scaling ranges can be formulated and solved, 
so that the solution tends to the exact one when the viscous and magnetic-diffusive cutoffs approach 
zero. In this approximation we are able to determine analytically the conditions for the existence of 
a dynamo effect and give an estimate of the dynamo growth rate. Among other things we show that 
in the large Prandtl number case the dynamo effect is always present. Our analytical estimates are in 
good agreement with previous numerical studies of the Kraichnan-Kasantzev dynamo by Vincenzi |14j . 

1 Introduction 

As we are lacking in the complete understanding of turbulence, it is useful to study simplified models which 
share many features of the full problem. One class of models are the passive advection models, where the 
velocity field is given some predetermined statistics. One then studies the effect of this velocity field on 
some other quantities such as the passive scalar, a temperature or a dye density in the fluid (see e.g. 0]). 
The term passive refers to the absence of backreaction to the velocity field. Naturally one would like to 
extend this study to passively advected vector fields. There exist physically realistic models such as the 
small resistivity magnetohydrodynamic equations for the magnetic field B G Mr interacting with a fluid 
(see e.g. El El d and references therein) 

<9 t B + v VB B Vv = kAB. (1.1) 
We will derive an equation for the pair correlation function 

(B i (t,T)B j (t,T')) (1.2) 

averaged over the velocity statistics, and attempt to solve it using a certain approximation scheme, which 
will be explained at the end of this introduction. The v above is the velocity field of a conducting fluid. It is 
assumed to be incompressible as in the constant density Navier-Stokes equation, that is, V-v = 0. Naturally 
the magnetic field is also incompressible due to the absence of magnetic charges, k is the resistivity divided 
by the vacuum permeability. The "smallness" of k here means that the equation is a good approximation 
when the magnetic Reynolds number Rm = VL/k is very large. Here L is the integral scale of the velocity 
field and V is the r.m.s velocity at such a scale. Since we assume from the beginning that L — > oo, the 
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approximation holds. In reality wc should also consider the backreaction of B on v, but we will only 
consider the passive case and just assume v to be given by the Kraichnan statistics (see e.g. 0] for a 
definition). It is defined as a Gaussian, mean zero velocity field with pair correlation function 



(v l (t,r)v J (t',r'))=5(t-t')D dk 



\k\ d +i 



/(L|k|)Py(k) =: S(t - t')Aj(r - r'; l v ) 



(1.3) 



with dk :— 



d d k 



and 



(1.4) 



to guarantee incompressibility. It is evident that Dij is homogenous and isotropic. The parameter £ for 
< £ < 2 describes the roughness of the velocity field with £ = 2/3 corresponding to the Kolmogorov scal- 
ing. The function / is an ultraviolet cutoff, which simulates the effects of viscosity. It decays faster than 
exponentially at large k, while /(0) = 1 and /'(0) = 0. For example we could choose f(l v k) = exp {—Ilk 2 ), 
although the explicit form of the function is not needed below. In the usual case without the cutoff function 
/ the velocity correlation function behaves as a constant plus a term oc r', but in this case we have an 
additional scaling range for r -C l v where it scales as oc r 2 . The length scale l v can be used to define 
a viscosity v or alternatively one can use k to define a length scale l K . We can then define the Prandtl 
number 1 measuring the relative effects of viscosity and diffusivity as P = v/k. Note that the integral scale 
was assumed to be infinite, i.e. there is no IR cutoff. 

In addition to an advection term v • VB in Eq. (|1.1|) . familiar from the passive scalar problem, there is 
also a stretching term B • Vv. This, and the fact that the magnetic field is a vector, will give rise to some 
interesting deviations from the passive scalar case. The most interesting one is probably the dynamo effect, 
an unbounded growth of the magnetic field depending on the roughness of the velocity field described by 
the parameter £, and the Prandtl number. This is in complete contrast to the passive scalar case, where in 
the absence of external forcing the dynamics is always dissipative 

The dynamo effect has been previously studied by e.g. Kazantsev where he derived a Schrodinger 
equation for the pair correlation function. However, the equation was still quite difficult to analyze except 
in some special cases. Some analytical and numerical results have been obtained e.g. in VS^ and ^3] 
(see the latter for further references). The goal of the present paper is to extend these considerations by 
introducing a set of approximate equations, which admit an exact solution. The analysis proceeds along 
the same lines as in a previous paper for a different problem by one of us [aj • The problem in the analysis 
can be traced to existence of length scales dividing the equation in different scaling ranges. In our case 
there are two such length scales, one arising from the diffusivity k and the other from the UV cutoff in the 
velocity correlation function. As will be seen in appendix A, what one actually needs in the analysis is the 
velocity structure function defined as 



This is all one needs to derive a partial differential equation for the pair correlation function of B, but it 
will still be very difficult to analyze. Hence the approximation, which proceeds as follows: 



1 We choose to write the Prandtl number as P instead of the usual Pr since it appears so frequently in formulae. 



X - ( ( Vi (t , r) - Vi (t, r')) (Vj (f , r) - Vj (f, r')) ) 




: 8{t-t')d l0 {T-v'-l v ) 



(1.5) 
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1) Consider the asymptotic cases where r is far from the length scales l K and l v with the separation 
of the length scales large as well. There are therefore three ranges where the equation is simplified into a 
much more manageable form. The equations are of the form dtH — M.H = 0, where M. is a second order 
differential operator with respect to the radial variable. We then consider the eigenvalue problem M.h = Xh. 

2) By a suitable choice of constant parameters in terms of the length scales, we can adjust the differ- 
ential equations to match in different regions as closely as possible. Solving the equations, we obtain two 
independent solutions in all ranges. 

3) We match the solutions by requiring continuity and differentiability at the scales /„ and l K . Also 
appropriate boundary conditions are applied. 

4) According to standard physical lore, the form of cutoffs do not affect the results when the cutoffs 
are removed. In addition to we can interpret l K as a cutoff. Therefore we conjecture that the solution 
approaches the exact one for small cutoffs. We also expect the qualitative results, such as the existence of 
the dynamo effect, to apply for finite cutoffs as well. 

For concreteness, suppose that M. is of the form 



The coefficients are some functions of the length scales l v and l K and the radial variable r. In general, 
solving the eigenvalue problem for such a differential equation is not possible except numerically. However, 
we can approximate the coefficients in the asymptotic regions when r is far from the length scales. The 
asymptotic coefficients are all power laws and solving the equations becomes much easier. Figure illus- 
trates this procedure corresponding to steps I) and 2) for any of the coefficients. 

After some preparations, we begin by writing down the equation for the pair correlation function of 
the magnetic field using the Ito formula. The derivation can be found in appendix A. The equation is of 
third order in the radial variable, but it can be manipulated into a second order equation by using the 
incompressibility condition. In section 2 the approximate equations will be derived when v <C k and n -C v, 
or Prandtl number small or large, respectively. We use adimensional variables for sake of convenience 
and clarity. The focus of the paper is mainly on the existence of the dynamo effect and its growth rate. 
Therefore we consider the spectrum of M . By a spectral mapping theorem, we relate the spectra of M. and 
the corresponding semigroup e* . It is then evident that if the spectrum of M. contains a positive part, 
there is exponential growth, i.e. a dynamo effect. 

1.1 Structure function asymptotics 

Due to the viscous scale l v in the structure function Eq. (|f .5|) , there are two extreme scaling ranges r 3> l v 
(incrtial range) and r <^l u . For r ^> l v we can set l u = in Eq. Ijf .5JI and obtain 



M = a(l„,l K ,r 



)d? + b(l„,l K ,r)d r + c(l u ,l K ,r). 



(1.6) 




(1.7) 



where 



r(i - £/2) 



(1.8) 



(d-l)(d + 2)' 



2 d+«-2 7r d/2 r ( d /2 + £/2)' 



The second case corresponds to the viscous range, which is to leading order in r: 
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Figure 1: A sketch of the procedure of approximating the example equation. The dashed vertical lines 
correspond to either one of the length scales l v and l K with pictures a) a plot of the "real" coefficient, which 
depends of the cutoff function (and is really unknown), b) an approximate form obtained by taking r far 
from the length scales (dotted parts of the lines are dropped) c) the approximations extended to cover all 
r e M and d) adjusting the coefficients to match at the scales l v and l K . For r much larger than the cutoffs, 
the error due to the approximation is lost. 



d<(v) := D 2 li- 2 r 2 ((d+l)fy - , (1.9) 

where 

We see that the viscous range form ()1.9JI can be obtained from ljl.7|l by a replacement £ — ► 2 and 
D\ — > Z?2^ -2 - Note that by adjusting the cutoff function / we can also adjust DijD\. 

1.2 Incompressibility condition 

Due to rotation and translation invariance, the equal-time correlation function of B must be of the form 

Gy(;,|x-x'|) := (B i (t,x)B j (t,x'))=G 1 (t,r)5 ij +G 2 (t,r)^-, (1.11) 
where r = |x — x'|. Additional simplification arises from the incompressibility condition diGij = 0: 

drGi - --^ I d r (r d - 1 G 2 ). (1.12) 
The general solution of the incompressibility condition can be written as 

Gi =rH'(r) + (d-l)H(r) 

G 2 — —rH'(r) ' 
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In terms of a so far arbitrary function H. Alternatively, adding the above equations we may write 



H(r) = ^- 1 (G 1 + G 2 ). (1.14) 

This observation leads to a considerable simplification in the differential equation for the correlation func- 
tion: whereas the equations for G\ and G2 are of third order in r, we can use the above result to obtain a 
second order equation for H . Then we would get back to G through Eqs. (| 1 . 1 H|> : for example we have for 
the trace of G: 

G u = (d - 1) (rH'{r) + dH{r)) , (1.15) 
although we refrain from doing this since H has the same spectral properties as Gu. 



2 Equations of motion 

The equation of motion for the pair correlation function is derived in appendix A: 

— Gij(t, r) = 2nAGij(t, r) + dapdj.ap — d a j^Gip^ a — dip^Gaj^ + dij^pGap. (2.1) 

The indices after commas are used to denote partial derivatives and we use the Einstein summation. By 
taking r 3> l v and r <C l v we can use the approximations l|1.7|) and (|1.9|) to write the equation in the 
corresponding ranges. This is done for the quantity H = (G\ + G2)/{d — 1) in the Appendix lAl as well, 
resulting in the equations 

^L=£(d- l)(d + £)D ir t- 2 H + [2(d + 1)k + (d 2 - 1 + 2£)£ir < n -H' 
at r 

+ [2k + (d - l)£>ir e ] H", r^>l v . (2.2) 

and 

^- = 2{d- l)(d + 2)D 2 li- 2 H + [2(d + 1)k + {d 2 + 3)D 2 lt- 2 r 2 ] -H' 
at r 

+ [2n+(d-l)D 2 li- 2 r 2 ]H", r «I„. (2.3) 
Simple dimensional analysis leads to the observation 



[k] = [IV £ ] = [D 2 lt 2 r% (2.4) 

where the brackets denote the scaling dimension of the quantities. We define the length scale l K as the scale 
below which the diffusive effects of n become important. This will be done explicitly below for different 
Prandtl number cases. In general, one can write k = Dilf,~ p lP for some p S (0,2]. Now one just needs 
to identify the dominant terms in the three scales divided by l v and l K . For sake of clarity, we choose 
to write these equations in adimensional variables. This can be done for example by defining r — Ip and 
t = 1 2 ~^t jD\ with I being a length scale. It turns out to be convenient to choose the larger of l K and l v 
as /. Since we deal with a stochastic velocity field with no intrinsic dynamics, we cannot, in principle, talk 
about viscosity. However, it is convenient to define a viscosity v (of dimension length squared divided by 
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time) by dimensional analysis from the length scale l v and the dimensional velocity magnitude D\ , giving a 
relationship between v, l v and D\ similar to what we would get in a dynamical model. We therefore define 



Dill 



(2.5) 



This permits us to define the Prandtl number in the standard manner as P — v / ml. We then consider the 
cases P < 1 and P > 1. 



2.1 Small Prandtl number 

I s- 



-M- 
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Figure 2: Sketch of the scaling ranges at small Prandtl number. 
Now v <§C At, and we choose as adimensional variables 

f r = l K p 
1 1 =3T T - 



(2.6) 



Note that the relation between l K and k has not yet been determined. In these variables the equations (I2.2|) 
and H2.3|) become 



d T H = f (d - l)(d + £)p- 2+ tH + 



2(d + l) 



(d 2 - 1 + 20/ 



D x l% 



+ (d-l)p i 



d 2 p H, p^l u /l R 



(2.7) 



and 



d r H = 2(d- l)(d + 2 



D 2 (I 



D x \l K 



£-2 



2(d+l)-^ + (rf 2 + 3) £>J 



D2 (lv\ 2 



d 2 p H, p<$Ll v /l K . 



Di \ l K 



5-2 



-d p H 



(2.8) 



As mentioned above, we also consider r -C Z K and r ^> l K , that is p <C 1 and p^S> 1, respectively. There are 
now three regions in p, divided by l v /l K and 1, with l v /l K <C 1. The regions, solutions and various other 
quantities will be labelled by S, M and L, corresponding to p <C l v /l K , l/i^^l and 1 -C p. See 
Fig. |21 for quick reference. Therefore the short range equation will be derived from Eq. <|2.8|l and the two 
others from (|2.7() . Consider for example explicitly the coefficients of d 2 p H: 



L 



M 



S 



K 

1 



+ (d-l)pt 



+ (d-l)pi 



(2.9) 



D 2 (lu\ _2 



! M + ,<i - 1) 5t lfl " 



G 



By definition of the length scale l K , in the region L the diffusivity is negligible and in the region M it is 
dominant, as it is in the region S since in there p approaches zero. The coefficients are then approximately 

L : {d-l)ffi 

M : (2.10) 

s - ; 



Dill 



Matching the coefficients of L, M at p = 1 provides us with a condition (matching between S and M gives 
nothing new) 

d-l=2-^-. (2.11) 

This is used as a definition of k as k = \{d — l)Dxl%. Writing down the short range equation with the 
above approximations, 

d T H s =2(d-l)(d+2)^(y\ H s + (d 2 -l)- p d p H s + (d^l)d 2 p H s , (2.12) 
by using the derived expression for the Prandtl number, 



P -H-7=T(V- (2 ' 13) 



and by defining 



D 



2 



D 1 Vd-1 



l-2/« 



(2.14) 



(remember that Z>2 could be adjusted by a choice of the cutoff function /, see Eq. ()1.9|l and below) a more 
neat expression is obtained for the short range equation. We can now write down all the equations: 

)±d p H s + (d-l)d 2 p H s (2.15a) 
\±d p H M + [d-l)d 2 p H M (2.15b) 
f 2£)pt- 1 d p H L + (d - l)p%H L . (2.15c) 

2.2 Large Prandtl number 

In/If 1 r> 

I § 1 M 1 L % 



d T H s = 


2(d- 


i)(d- 


-2)P 




f (d 2 


9 t Hm = 


e(d- 




-Op" 


- 2+ ^H M - 


h(d 2 


d T H L = 






-Op" 


2+ ZH L + 


(d 2 - 



Figure 3: Sketch of the scaling ranges at large Prandtl number. 
Now 2/ 3> /c, and we choose 
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Then the equations l|2.2[l and (|2.3(l for r ^> l v and r <^l v become in the new variables 



(2.16) 



d T H = t(d-l){d + (,)p- 2+ tH 



2(d+l) 



d 2 p H, p»l 



(d 2 -l + 2^ 



-d p H 



(2.17) 



and 



d T H = 2(d — l)(d + 2)—H + 
D\ 

k , , s Do o 
' 2 ^ + {d - 1] D-/ 



2(d + l) 



Dili 



d 2 H, p«l. 



2 „2 



(2.18) 



The ranges S, M and L now correspond to p <C I k J lv «p< 1 and 1 -C p, see Fig- EI Note that 

equations in both S and M are now derived from Eq. I|2.18[) . As before, we consider again the coefficients 
of d 2 H and drop the terms oc k in L and cx p 2 in 5. The diffusive effects are not dominant in the region 
M since r 3> l K , so we drop the oc k term in M too. The approximative coefficients are then 



L : (d- l)p e 
M: (d-l)^p 2 

s 



Dili 



(2.19) 
(2.20) 



We then obtain two equations by matching the coefficient of L with M at p = 1 and of M with S at l K /l v : 



(d-1) 



£>2 

D 2 



(d-l) = (d-l), 

2 



= 2- 



(2.21) 



with solutions 



D 2 = D l 

K = ^D.fjt 2 . (2.22) 

The Prandtl number is in this case 

2 (l^ 2 
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d T H s = 


2(d- 


l)(d- 


-2)H S + 2^±d p H s + ld 2 p H s 


9 t Hm = 


2(d- 


l)(d- 


- 2)H M + (d 2 + 3)pd p H M + (d- l)p 2 d 2 H M 


d T H L = 






- t)p- 2+ ZH L + (d 2 -l + 2^- 1 d P H L + (d 



Note that one can obtain this from the small Prandtl number equation l|2.13[) by replacing £ — > 2. This is a 
reflection of a more subtle observation that the large Prandtl number case for any £ is similar to the small 
Prandtl number case with £ = 2. We collect the equations using the above approximations, 

(2.24a) 
(2.24b) 

)ptd 2 H L . (2.24c) 

Note that the short and long range equations are somewhat similar to the respective small Prandtl number 
ones, Eqs. H2.15al2.15bl2.15cjl . However, the equation in the medium range above is scale invariant in p, 
unlike the corresponding small Prandtl number one. 

3 Resolvent 

Given a differential operator M. with a domain D (A4), we define the resolvent 

R(z,M) := {z-MY 1 (3.1) 

and the resolvent set as 

p(M) := {z e C\z - M : D{M) -> X is bijective} . (3.2) 

The complement of the resolvent set, denoted by a(A4), is the spectrum of M.. According to the Hillc- 
Yosida theorems (see e.g. [2]), if (A4, D(M)) is closed and densely defined and if there exists zo€K such 
that for each z€C with 5Rz > zq we have z £ p(M), and 

\\R(z,M)\\<-^—, (3.3) 

MZ — Zq 

then M. is the generator of a strongly continuous semigroup T(t) satisfying 

|| T(t) || < e Zot (3.4) 

and vice versa. If in addition M. is a so called sectorial operator, meaning that its spectrum is contained in 
some angular sector {z £ C : | arg(z — z )| > a > ir/2} and that outside this sector the resolvent satisfies 
the (stronger) estimate 

\\R(z,M)\\<-^-, (3.5) 
V- — zo\ 

then M. generates an analytical semigroup, for which the spectral mapping theorem 

a(T{t)) = {0}Ue ta{M K (3.6) 

holds, relating the spectrum of the generator to that of the semigroup. One can also use the Cauchy integral 
formula 
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T(t) := e tM = — [ dze zt R(z,M), (3.7) 
2m J 

c 

where the contour surrounds the spectrum a(M). 

We do not prove here that M is sectorial, however we refer the interested reader to the general math- 
ematical theory in where it is explained and substantiated that strongly elliptic operators are, under 
quite general assumption, sectorial generators, on a wide range of Banach spaces (e.g. LP and C 1 spaces to 
name but a few). 

According to the above discussion, in order to explain the existence of the dynamo effect and its growth 
rate, we only need to find the spectrum of M via the resolvent set p{M). Note that we are interested only 
in the positive part of the spectrum, since we want to determine the existence of the dynamo effect only. 

The operator M in our case is cut up as the operators Ml, Mm and Ms in the corresponding ranges, 
obtained from the equations (|2.15al l2~15bl I2.15cfl and 12.24al l2~24bl I2.24cf) . The resolvent is found from 
the equation 



(z-M)R(z,M)(p,p') = S(p,p'). (3.8) 

Since we are primarily interested in the long range (L) behavior p > 1, we let p' stay in the region L at all 
times. This results in three equations 

(z-M L )R L (p,p') =8{p-p') 

(z-M m )Rm(p,p') =0 (3.9) 
(z-M s )Rs(p,p') =0 

where Rl(p, p') is the expression of the resolvent for p S L (the large scale range) and p 1 G M + and similarly 
Rm and Rs are valid when p is in the middle and small scale ranges respectively. We require the following 
boundary conditions from the resolvents: for small p we are in the diffusion dominated range, so we require 
smooth behavior at p — > 0. For large p we eventually cross the integral scale (although we haven't defined 
it explicitly) above which the velocity field behaves like the £ = Kraichnan model, leading to diffusive 
behavior at the largest scales for which the appropriate condition on the resolvent is exponential decay at 
infinity. 

3.1 Solutions 

Assuming p ^ p', we solve the equations l|3.9|l with the corrsponding operators M. 

The operator Ml does not depend on the Prandtl number. So in the region L, we get from e.g. 
Eq. H2.24cfl (we use lowercase letters h^ 1 to denote the independent solutions) 

ht( P )= p - d/2 -^Z±(w P ^), (3.10) 

where Z^ = I\ and = K\ are modified Bessel functions of the first and second kind respectively, and 
we have introduced 



(d-l) 

and the order parameter A is 



(3.11) 



_ y/d[2{d-l)*-(d-2){2ti + d=m 
X (2-fl(d-l) • (3 - 12) 



10 



Because the range S is always in the diffusive region, we require smoothness of the solution at zero. 
Only one of the solutions satisfies this, so we get from Eqs. H2.15al) and H2.24a|) 



h SA (p) = P~ d ' 2 h/2 ( Jj^j - 2(d + 2)pi-2/£p ) (3.13) 



and 



HsAp) = P~ d/2 Id/2 (J 7 --{d-l){d + 2)^P p) , (3.14) 

where the subindex 1 refers to P <C 1 (small Prandtl number) and 2 to P ^> 1. We will use this notation 
in other objects as well. 

In the range M, when P> 1 we have the scale invariant equation in (|2.24b|) with power law solutions 

hi I>2 {p)= P - d/2 -^ ±i , (3.15) 

where 



C - ^ (3.16) 

with 

z 2 = -^[(2 - £)A]\ =2 = _ i ^-^(d 3 ~ 10d 2 + 9d + 16) . (3.17) 

The medium range equation for P <C 1 cannot be solved exactly, but we can consider it in two different 
asymptotic cases. From Eq. (I2.15b|) we get 



£(<* + £V~ 2+? - Rm + {d+ l)- p d p R M + d 2 p R M = 



(3.18) 



and note that since by definition of the medium range l v jl K <p<l, i.e. 1 < p~ 2+ ^ < {l K /lv) 2 ~^ (the -C 

was replaced by <, so that things remain valid even as £ — > 2), the term oc p~ 2+ ^ can be dropped if we 
assume that 

\z\ > (WM 2 -« « . (3.19) 

If on the other hand we have 

\z\ <1, (3.20) 

then z can be neglected in the equation. 

The solution for large z is similar to the short range solutions, 



= p- d,2 Z% 2 (Jj^p) » 1*1 > (lM 2 ~ e , (3-21) 
where we denoted the P -C 1 case by a subscript 1. For small z we have instead 

h% = p- d/2 Z^ (2Vd7? + V /2 ) , 1*1 < I- (3.22) 

where Z^/i = *^d/£ an< ^ ^ are Vessel functions of the first and second kind respectively. It turns 
out however that the explicit form of the above solutions affects only a specific numerical multiplier and 
has no effect on the presence of the dynamo. Because of this we in fact derive a lower bound for the growth 
rate which in view of the present approximation provides a more reliable result. 
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3.2 Matching of solutions 



Consider equations l|3.9[l . We denote the long range regions p < p' and p > p' as L < and L>. The boundary 
conditions for the resolvent demanded finiteness at p = 0, but in general the resolvent must be in £2(1$+). 
We therefore have in the region L> only the solution, since it decays as a stretched exponential at 
infinity (the other one grows as a stretched exponential) . We also drop the subscripts labelling the different 
Prandtl number cases for now. The full solutions are written as follows: 



Rs{z\P-P) =ah s (p) 

(3.23) 



R M {z\p, p') = C^h+jip) + C M h M (p) 



RlMPiP') =cth + L {p) + czh- L { P ) 
R L> (z\p,p') =Phl(p), 

We denote the matching point between the short and medium ranges by a.;, i.e. 



{:=t/t <" 4 » 

The other matching points are p — 1 and p = p' in both cases. There are six coefficients to be determined, 
a, C^, and (3, and in total six conditions, four from the continuity and differentiability at p — di and 
p = 1 and two conditions at p = p' around the delta function, so all coefficients will be determined from 
these. They will then depend on the variables z and p' . The C 1 conditions at p = 1 are 

C+h+(l) + ClhlO) = C+h+(l) + C M h M {\) (3.25) 

and 

C+dh+(l) + C- L dh- L (l) - C+dh+(l) + C M dh M (l), (3.26) 
where we denoted dh(l) — d p h{p)\ p=1 . This can be expressed conveniently as 



ht h- L \ fC+\_fh+ M h M \ fC+ 



dh+ dh L J 1 \C L J \dh+ dh M J \C M 



(3.27) 



where the matrix subindex refers to evaluation of the matrix elements at p = 1. Since we have only one 
solution at short range, we get similarly at cti 



h + M h M \ (C+\ =a (h s 
dh+j dh M ) \C M J \dh s 



M =a "°J • (3-28) 



where again the matrix subindex indicates the point where matrix elements are to be evaluated. We can 
solve these for C^, 



\ C lJ \~ dh L h L J 1 \ dh M 9h Mj 1 \- dh M K J ai \ dhs J a t 

The numeric constant J' above contains the determinants of the inverted matrices and a. It is certainly 
nonsingular due to the linear independence of the solutions. We have decided not to explicitly write it 
down since, as we will see below, we only need the fraction C^/C^. Now we have piecewise the resolvents 
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( r l< ( Z \p,p') =ci(hi{p)+%?hz(pj) (3 30) 

\ R L> (z\p,f/) =0hi(j>), 
and we still need to use the first equation of l|3.9[l for and 0. The continuity condition is 



C+[ht(z,p') + ^hl(z,p')) =phj,(z,p'). Cl.TJ.) 



The other condition is obtained by integrating the equation with respect to p over a small interval and then 
shrinking the interval to zero: 

C+ (dh+{p') + ^dh- L (p'^j - Pdh- L {p') = 1. (3.32) 

These can be solved to yield 



C + = -h¥±- - (3.33) 

w(hi,hz)( P >) 

and 



C+ht{pi) + Clh- L {pi) 

13 - ctw(ht,h- L )( P >) ' (3 - 34) 

where W is the Wronskian, W(f,g) = fg' — fg. Explicitly from Eq. I|3.1()|l . 

W(ht,hl)(z,p') = (p')- d -^W(h,K x ) = -{l-Z/2){p')- d - l -A. (3.35) 
Using the above obtained expressions of C~[, (3 and W in Eq. H3.3()(l we thus have the solutions 



RlMPiP') = - LE^^Zll (h+{p)h- L {pi) + (§0 h- L { P )h- L { P ')) 
R L> (Z\P,P') = - (p ' ) T-€/2 /( '" 1) (hKp)hi(p') + (g) hl{p)hlW)) 



(3.36) 



with C L /C^ obtained from Eq. (|3.29|) . We have calculated in appendix B the asymptotic expression for 
C^/C^ for the two Prandtl number cases P <C 1 and P ^> 1, 

= dh+(l) - A L h+(l) 
C+ dhl(l)-A L hl(l) [ ' ' 

with leading order contribution to as 



dhtr(l) 

n M\ l ) 



with cither the + or the — solution understood, the choice depending on the Prandtl number. 
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4 Dynamo effect 



The mean field dynamo effect for the 2-point function of the magnetic field corresponds to the case when 
the evolution operator M. has positive (possibly generalized) eigenvalues. Eigenvalues correspond to poles, 
in z, of the resolvent given in Eq. H3.36[l and generalized eigenvalues to branch cuts. The z dependence 
is not seen explicitly in Eq. (|3.36|1 . but recall from Eq. 1)3. 10[) that the hj^ depend on z, and we see from 
Eq. H3.37[l and matter in Appendix iBl that there is further dependence through /C^ . 

The hf; (see Eq. (|3.10fl % ) depend on the square root of z, and since the Bessel functions are analytic on 
the complex right half-plane, this square root dependence leads directly to a branch cut along the negative 
real axis in the z dependence of the resolvent. This corresponds to a heat equation like continuum spectrum 
of decaying modes, these modes don't contribute to the dynamo effect. 

Any other possible contributions to the spectrum come from the fraction C^/C^. An expression for 
the latter is given in Eq. (|3.37|) . with A^ computed in Appendix IBl Equation (|3.37[1 can be simplified by 
noting that (using Eqs. (|3.10|) and l|3.11|l l 



dhfil) = - (I + ^y) ZfM + (1 - £/2Ri±« (4.1) 

Then we can write 



CI 



(1 


- £/2) W I' x (w) - 


2 




h(w) 


(1- 


Z/2)wK' x {w) - 


d _ 
2 


- A+ A ^_ 


K x (w) 



(4.2) 



where we have replaced the partial derivative symbols with primes. We underline again that w depends on 
the square root of z, so the complex plane minus the negative real line for z corresponds to the ^Rw > 
half-plane for w. Since Bessel functions are analytical on this half-plane, the new singularities introduced 
by C^/C^ may come either from singularities of A^ or zeros of the denominator . The latter (i.e. 
= 0) may be written more conveniently for future needs as 



w 



K' x (w) _ 2 
K x (w) ~ 2~e 



l 



A, 



(4.3) 



where we define A^. 

Finally we remind the reader that z corresponds to the growth rate with respect to the reduced time r 
rather than real time t, and the real growth rate is thus (r/t)z, where for the P — > case from Eq. (|2.6(l 
we have r/t = Dil\~ 2 and for the P — > oo case from Eq. I|2.16|l we have r/t = D\l\^ 2 . 



4.1 Prandtl number P —> 

Based on asymptotics in Appendix lB.il we expect that A^ does not introduce new singularities in the small 
Prandtl case, so we only need to deal with Eq. (|4.3(l . 

Let us first study the large w asymptotics of the two sides of Eq. (|4.3(l . As shown in Appendix IB. II 
Eq. (|B.18J) . in the limit of vanishing Prandtl number and for large z, - or equivalently large w- from the 
medium range solutions of Eq. (|3.21|) we get 



id/2 ((i - CAM 

For large w we deduce from the asymptotic properties of Bessel functions [Zj that A^ ~ (1 — £/2)w. In the 
same manner we deduce that the left hand side of Eq. Ij4.3|l is wK' x (w) / K\(w) ~ — w. 
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Next we remark that if A is pure imaginary, which based on the definition of A, Eq. 1)3. 12)1 . happens for 

e>r:=c^-D ? (4.5) 

then K\(w) has an infinity of positive zeros (may be seen from its small w development), accumulating at 
w = 0, and more importantly it has a largest zero (may be seen from its large w asymptotics) , which we 
shall denote by wq. At each zero of K\(w), the l.h.s. of Eq. 1)4. 3)1 has a pole. 

Since for large ui, asymptotically K\(w) ~ (-k/2w) 1 ^ 2 exp(— w) > 0, we must have K' x (wq) > (note 
that we can exclude K\(wq) = K' x (wq) = 0, since the Bessel function is the solution of a second order 
homogeneous differential equation). Hence the pole of w K' x (w) / K \(w) at wq has a positive coefficient. 

From the above asymptotic properties of the two sides of Eq. 1)4.3)1 . in conjunction with the positivity 
of the coefficient of the pole at wq, using continuity of the two sides we deduce - see also Fig. 0]- that 
Eq. 1)4.3)1 admits a solution which is at least as large as the largest zero wo of K\(w). 

As a convenient means of estimating the solution of Eq. 1)4. 3)1 . we may use the lower bound wq. One 
should also be able to obtain an upper estimate, which we expect to be of the same order of magnitude as 
wo, as argued below. For the d = 3 dimensional case we plot wq in Fig. [3] and it indeed compares well with 
the numerical results of |14) . the latter involving no approximations. This corroborates our approach. 

First we note that, as shown in Appendix lC.il A^ is increasing. The small z asymptotics of Al is given 
in Eq. 1)4.6)1 . Two cases are distinguished, depending on wether Al(z = 0) > or not. If it is then we 
have the upper bound w' , where w' Q is the largest zero of K' x . Otherwise we use (wK' x (w) / K \(w))' < — 1 
from Appendix IC . 21 and get the upper bound w\ — w' Q + \ Al(z = 0)|. Fortunately for £ near £* we have 
Al(z = 0) > (can be seen by numerical computation for relevant values of d), which means that the 
£ — > £* asymptotics computed in Sect. 14.21 is valid. The situation for the £ — * 2 asymptotics is more 
complicated, at the end of Appendix IB . 1 1 it is explained why near £ = 2 we may use the upper bound w' : 
although Al(z = 0) < 0, we can justify Al(z) > for z the dynamo growth rate. 

We can also study the small w asymptotics of the two sides of Eq. ()4.3)) . As shown in Appendix IB. II 
Eq. I)B.21)I . in the limit of vanishing Prandtl number and for small z, - or equivalently small w- from the 
medium range solutions of Eq. 1)3.22)1 we get 



-ji/ C +i(2Vd/e+i) 



Al = -g A /d/m , ■ ( 4 - 6 ) 

J d/i {2^dW+l) 

which is obviously independent of w, so it is in fact the w — > limit of A^. Finally we deduce from the 
asymptotic properties of Bessel functions jjj that, as w — > 0, if A is real then the left hand side of Eq. 1)4.3)1 
goes to —A, and if A is pure imaginary then it has an infinite number of poles at the positive zeros of K\. 

So far we have dealt with the case when £,*<£,< 2, and we were able to show the existence of a dynamo 
effect and give a lower bound on the dynamo growth rate. We are now going to argue that in all other 
situations there is no dynamo effect. 

Using the integral representation of the modified Bessel function of the second kind, 



oc 

K\(w) = Jdte- wcoah ^cosb(Xt), (4.7) 
o 

we see that when A 6 I (equivalently when < £ < £*) and w > 0, then K\(w) > 0. On the other 
hand using the recurrence relation K' x (w) — —K\-i(w) — XK\(w)/w we get wK' x (w)/K\(w) = —A — 
wK\_i{w) / K\{w) < —A, where the last inequality follows from the positivity of K\-i(w) and K\{w). 
Using the power series development of Bessel functions we also get that wK' x (w)/K\(w) — > — A as w — > 0. 
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Figure 4: Based on the asymptotic properties of the left and right hand sides of Eq. 1(4.3(1 . when A is pure 
imaginary and K\(w) has a largest zero wq, the latter has to be a lower bound on the largest solution of 
Eq. 1(4.3(1 . The dashed line depicts the right hand side of the equation. 

For 3 < d < 8 the critical £* takes values between 1 and 2, in particular for d = 3 we get £* = 1 as 
expected ^lElE]- In this case for £ such that < £ < £* < 2, the right hand side of Eq. 1(4.3(1 is always 
positive (can be seen numerically/graphically), so there can be no solutions of Eq. (|4.3(l . hence no dynamo. 

For d = 2 one readily verifies that £* = +oo so there is no "normal" dynamo (i.e. one with £* < £ < 2), 
and one further checks that there isn't any "exceptional" solution with £ < £* (cf. Sect. 15.3(1 . by showing 
that the left hand side of 1(4.3(1 starts from — A at w = and is decreasing, while the right hand side starts 
from a larger value and is increasing. 

Finally, for d > 9 we get £* > 2 so there is no "normal" dynamo either, and plots indicate that the left 
hand side of Eq. (|4.3|) is always negative while its right hand side is always positive, so again there is no 
dynamo at all. 

We have thus found that in dimensions 3 < d < 8 a critical value 1 < £* < 2 exists above which the 
dynamo is present and below which we don't expect it to be present. In other dimensions we expect no 
dynamo for any value of £. 

4.2 Asymptotics for £ near £* and 2 

In this subsection we give estimates of the growth rate of the dynamo in the cases when £ is near the 
critical value £* above which dynamo is present, and when £ is near its maximum possible value 2. What 
we need is an estimate of the largest solution w of Eq. (|4.3|l . from which the corresponding growth rate z is 
immediately deduced through Eq. 1(3.11(1 . As an order of magnitude estimate for the solution, the simplest 
is to take the largest zero w of K\{w), as can be seen by inspecting Figure 0] 

The case £ \ £*, corresponding to A — > along the imaginary axis, is somewhat simpler, it can be 
dealt with starting from the integral representation 1(4.7(1 . Since £ > £*, the parameter A is imaginary 
and we write A = iX with A G R, hence K^w) = J dt exp(— w cosh(t)) cos(Ai). Now cos(Ai) is positive 
near t = and it becomes negative for the first time only for t > tt/(2X). On the other hand the term 
exp(— w cosh(t)) is basically a double exponential and decays very fast for wcosh(t) > 1. So in order to get 
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Figure 5: In the above figure we have plotted the lower bound wq with the numerical results of |14j . The 
middle figure shows plots of both data with 1/ log(e(£)) 2 on the y-axis. The lower bound data shows linear 
behavior consistent with the asymptotics in eq. 1)4. 8[l . We note that there seems to be a numerical error 
in the data of for £ = 1.02. In the lowest figure we have also plotted the numerical data in near 
£ = 2 for (15/2 — e(£)) 3 ^ 2 showing linear behavior as expected in eq. 14.91 The plots in other dimensions 
look similar, except that they begin from the critical value £* > 1. 
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for the previous integral a non-positive result, we need wq cosh(-7r/(2A)) ~ 1 implying wq ~ exp(— 7r/(2A)). 
Through Eq. (|3 . 1 1 1> one deduces the behaviour Inz ~ c/A, and since A oc - S,*) 1 ^ 2 near £* (the term 
under the square root in Eq. (|3.12() is expected to have a simple root at £ = £*), we finally have 

lnzoc (£-fT 1/2 (4.8) 

near £*. 

We now pass to the asymptotics of the case £ — > 2. Under this limit A diverges as (2 — along 
the complex axis. From plots or asymptotical formulae we can convience ourselves that the largest zero of 
K\(w) occurs at w ~ |A|, so that is the region where we are going to look for it. We write K in terms 
of the Hankel function of first kind K x {w) = fe mX / 2 H^\iw) and use the approximation formula in the 
transitory region (i.e. when parameter and argument of the Bessel function are of same order): 

H x (A + tX ) ~ 3A i/3 t H i/3(— t ) ■ 

We deduce that there is some to > such that w — | A | — i 1 A | ^^Z 3 = |A|(l-i |A|~ 2 / 3 ). Combining Eqs. H3.11j) 
and 13.12JI we get 

*ra 2 2 - c(2 - £) 2/3 , (4-9) 
valid for £ near 2, where c > is some constant of order unity and Z2 was introduced in Eq. (|3.17() . 



4.3 Prandtl number P — > oo 

For large Prandtl number the analysis proceeds exactly as in the previous section, except that we have a 
different A/,. From Eqs. (|B.30(I and (|B.27|I we have, using the definition of £ from Eq. (|3.16|l . 



**-<-jpt-5- <"°> 

There is now a branch cut originating from Z2 (defined in Eq. (|3.17|l ^ extending to infinity along the negative 
real axis. When 3 < d < 8, Z2 is positive and the branch cut extends up to the positive value Z2 along 
the positive real axis, i.e. the spectrum has a continuous positive part for all £ £ [0,2]. Another major 
difference when comparing to the small Prandtl number case is that the spectrum is continuous also in the 
positive part. 

We conclude that the dynamo is present for all £ for large Prandtl numbers. 



5 Some remarks 

5.1 Schrodinger operator formalism 

We have performed our analysis in the diffusion process setting, but it could have equally well been done 
in the Schrodinger operator formalism (e.q. |14|). basically with the same kind of cutting up and piecing 
together technique. 

Consider the zero energy Schrodinger equation 



iP"(r) = V(r)ip(r), (5.1) 

where V(r) = m(r)U(r) is the effective potential. The potential V behaves as 2/r 2 at very short and 
long scales, but as — 4/r 2 at the medium range. The medium range solutions are sin (VT5/21og(r)) and 
cos(vT5/2 log(r)). When the Prandtl number is increased, the medium range region is stretched, and it is 
clear that for sufficiently large Prandtl numbers the solutions cross zero an increasing number of times (see 
Fig.[fjJ). According to a well known theorem, such a solution cannot be a ground state (see e.g. PP- 90). 
In fact, the number of zeros of the solution (with nonzero derivative and excluding the zero at r = 0) is the 
number of negative energy states, which implies the existence of unbounded growth. 
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Figure 6: Sketchy plot of the medium range zero energy solution sin(vT5/2 log(r)) crossing zero, implying 
the existence of a negative energy state. The dashed lines correspond to regions where the middle range 
approximation is no longer valid. As the range grows when increasing the Prandtl number, more zeros will 
emerge. 

5.2 Finite magnetic Reynolds number effects 

Let us finally touch upon some questions not discussed in the text. Our method allows us in principle, 
without further complications, to estimate the critical magnetic Reynolds number (dependent on velocity 
roughness exponent £ and space dimension d) at which dynamo effect sets in, and the growth of the dynamo 
exponent with Reynolds number. However we get only a logarithmic estimate whose uncertainty is at least 
an order of magnitude or even two, which makes it not too useful. Notwithstanding, we would like to 
mention that the estimates we would obtain this way are hardly compatible with numerical results of |14j . 
our thresholds being significantly lower. This issue is currently clarified with D. Vincenzi. 

5.3 Exceptional solutions 

An other issue is that of the existence of "exceptional" dynamos. It seems to us that the "typical" dynamo 
(note that we consider here only the infinte magnetic Reynolds number case) corresponds to the situation 
when our £ > £*, in which case there is an infinite discrete spectrum of growing modes. However our 
equations do not exclude a priori the possibility of a single growing mode at some £ < £*. In fact, if 
we take for example, at a formal level, d = 2.125 then £* « 1.82 and for (where £' is some 

value of which we only need to know here that £' < 1.77), Eq. I|4.3[l will have, in what we have called the 
small z approximation (cf. Eq. 1)3. 20)1 ). a single solution wo > 0. If we take £ = 1.77 then wq « 0.077 and 
z w 8.8 • 10 < 1 in a self-consistent manner. However it remains to be known if such a solution is not 
just an artefact of our resolution method, and if not, then to see if one can construct a model where such 
solutions occur for the more physical value of d = 3. 

6 Conclusions 

The mean-field dynamo problem was considered in arbitrary space dimensions. We have shown that, to 
obtain the spectrum of the dynamo problem, the equation (|4.3[) has to be solved for w, from which the 
growth rate z can be expressed through Eq. (|3.11fl . The quantity A^ appearing in Eq. I|4.3|l is given, for 
small magnetic Prandtl numbers, by either Eq. 14.4JI or 1)4.6(1 . depending on which of the self-consitent 
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conditions l|3.19[) or l|3.2U[l is verified (note that this leaves a gap between with no formula). For large 
magnetic Prandtl number we have to use Eq. I|4.10|l instead. 

It was observed that, in our model, the dynamo can only exist when 3 < d < 8. The results for small 
Prandtl numbers were shown to confirm previous results (141 l§] obtained in three dimensions. For d > 3 a 
critical value for £ was found, above which the dynamo is present, which is larger than the three dimensional 
critical value £* = 1. Furthermore, in the vanishing Prandlt number limit we have obtained the asymptotic 
estimates (|4.8|l and (|4.9fl . which are in good qualitative agreement with numerical simulations of |14| . 

For large Prandtl numbers it was shown that the dynamo exists for all £ and that the spectrum is 
continuous. We hope our work will contribute to clarifying this somewhat controversial issue. The physical 
idea behind our explanations is that at large magnetic Prandtl number the magnetic field can feel the 
smooth scales of the fluid flow (they are not "wiped out" by magnetic diffusivity) , and correlations in the 
velocity field above the magnetic diffusive scale l K won't do more harm to the dynamo than if we had a 
Batchelor type flow with no correlations of velocity at scales significantly larger than l K . 

Our methods were based on approximating piecewise the evolution operator of the two point function 
of the magnetic field. This approximation introduces inaccuracies and one may ask how these influences 
the fine details of our reasoning, which relied on not so evident estimates. We think that the general 
picture sketched up should be valid for the exact problem also, based on the good agreement with available 
numerical data from the litterature. Since for £ = £* and £ = 2 one can find the fastest growing mode 
explicitly, it should also be possible to do a perturbation theory around these points and place our results 
on a firmer ground. This is however left for future work. 
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A PDE for the 2-point function of B 

We rewrite equation Ijl.lfl as an Ito type SPDE (following the formalism of ^UJj sec t- 5) 



dBi + dw ■ VB, - B • Vdw l - k' AB.dt = 0, (A.l) 
where k! = k + D/2 with D defined as 

D ij (Q)=DS ij . (A.2) 

The new diffusion term in k! emerges by advecting the magnetic field along the particle trajectories similarly 
as in the passive scalar case by using the Ito formula. It will cancel out eventually, as it should. We can 
express the above equation more conveniently by defining 



dh(t, x) = -V x ijk {Bj(t, x)dw k {t, a?)) , (A.3) 

_ X.. Fix 2 rp U ^ 

ijk z 3 k 



where T>f ik = Sijd^ — Sikd* ? The equation is then simply 



dB t ~ n'AB.dt = dbi. (A A) 

For a function F of fields B, we have the (generalized) Ito formula, 

f 5F If 5^F 

dF(B(t,.))= dx——[ K 'AB l dt + db l ] + - dxdy - E (db i (t,x)db j (t,y)) . (A.5) 

J 6Bi(x) 2 J dB i (x)dB j {y) 

The advecting velocity field is a time derivative of a Brownian motion on some state space, that is 

~Edwi(t,x)dwj(t,y) — dtDyix — y), (A. 6) 

where was defined in Eq. This means that 

Edh(t, x)db j (t, y) = V x m V v jmn (u k (t, x)u m (t, y)D ln {x - y)) . (A.7) 

We apply this to F = Ui(t,x)uj(t,y), denote Gij(x — y) = ~Eiii(t, x)iij(t,y) and use the decomposition 
Dij(x — y) = DSij — dij(x — y) introduced in Eq. Noting that terms proportional to dw disappear, 

we obtain the equation for the two point function: 

—Ga (t, r) = 2nAGij (t, r) + , (A.8) 

where r = x — y and 

J~ij — d a ^Gij^ a p d a j^Gi{3^ a di{3 j0 ,G a j^ 

f) + di^apGajj (A.9) 

where the indices after commas denote partial derivatives. Note that this depends only on k, not k', i.e. 
the constant part DS^ of the structure function is absent. Using the decomposition (|1 . 1 1ft and the explicit 

2 This is just a rewriting of the expression V X (B X v) for incompressible fields B and v 
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form of the long distance velocity structure function (|1.7I) wc get from Eq. I|A.8|) two equations for Gi and 
G2, 



(d-l)rG[ + r 2 G'l)+A, 



(A.10) 



and 



^ = (_2dG 2 + (d - l)rG' 2 + r 2 G'i) + 6. (A.ll) 
at r z 

The symbols ^4 and B are the terms arising from the interaction with the (long distance) velocity fields. 
Using the relations (jl . for G\ and G2 in terms of H, their explicit form is as follows: 

A £(d-l)(d-3 + 0ff(r) 



D ir -2+C 

. . . + (2 - d - 2d 2 + d 3 + (-5 + d + 2d 2 )£ + (1 + d)£ 2 )rH'(r) 

... + (2d(d - 1) + (d + l)0)r 2 ff"(r) + (d - l)r 3 ff"'(r), (A.12) 



B 



= -£(d-l)(2-0ff(r) 

+ ((1 - d 2 ) + (d - 5 + 2d 2 )£ + 4£ 2 - £ 3 )rH'(r) 
+ (d+ l)(d - 1 + t)r 2 H"{r) + (d - l)r 3 H"'(r). 



(A.13) 



Now we can just add the equations l|A.10Jl and 1|A.11J| . and by using Gi + G2 = (d — we get 



^ = ^(d-l)(d + 0Ar- 2+ «ff(r) 
. . . + (2(d + 1)« + (d 2 - 1 + 2f) J D 1 r ? ) r^H'ir) 

... + (2K+(d-l)Dxr £ )H"{r). (A.14) 
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B Computation of the fraction C L /C£ 

By evaluating the matrix multiplications on the right hand side of Eq. 
C7/C+ as 



we can write the fraction 



CI = dh+{\) - A L h+{\) 

c+ dh- L {\) - K L h- L {iy 

where can be written as the following nested expression: 



(B.l) 



Am 
As 



dh+ (l)+A AJ 3fe- (1) 
h+ (l)+A A fh M (l) 

(q,)-A s h+ (a,) 
9/t^(ai)-As/i^(o») 

hs{a,i) 



(B.2) 



This follows from defining 



h S 
dh S 



hs(a.i) 



A, 



(B.3) 



and writing equation 



as 



dh L —h L 
-dht hf 



h + h~ 

dh+j dh~ M 



dh iri 

-f)h + h + 



M 



'M 



A s 



(B.4) 



where hs(cti) is absorbed in the coefficient. We have defined above a constant c which gets cancelled in 
the end of computations. It will be used below as well as a generic constant that does not affect the final 
results. Multiplying the last matrix with the vector, we define similarly 



dh M h M 

un M IL M 



I 

As 



dh M - A s h M 



-dh 



M 



Ash 



S a M 



1 

A M 



(B.5) 



that is, 



A M = - 



dhU 



Ash^ai) 



dh M (di) - A s h M {a,j) 
Doing this again for the second matrix, we obtain similarly 



(B.6) 



A, 



and finally 



dh+(i) + A M dh M (i) 

/i+(l) + A M /i M (l) 



(B.7) 



CI = dh+(\) - A L h+{\) 
C+~ dh£(l)-A L hi(lY 
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We are interested in the leading order behavior of the fraction C£ / only, so we need to determine 
what happens to Am as P approaches zero or infinity. It turns out that either Am -> or Am -> ±00, so 
to the leading order, 



B.l P<1 



Below the suspension dots denote higher order terms in powers of P (or P 1 for large Prandtl numbers). 
Recall from Eqs. jjQUl and [(27131) that 

d-1 ^ 



ai = l v /l K = (-^- p ) ■ ( B - 10 ) 
The short range solution was 

hs(p) = P~ d/2 h/2{a P ) (B.ll) 



with a temporary notation a = \J [z + 2P 1_2 /?(2 — d — d 2 ))/{d — 1) and note that |a| behaves as P 1 / 2 
Using standard relations of Bessel functions |Zj and using the definition for A5 in Eq. (|B.2p . we have 

As — — — p — - — a t t-t- — r— . (B.12J 

^ ai ) I d/2 ((^P) m a) 

Since P is small and the arguments of the Bessel functions above scale as P 1 / 2 , we can use the expansion 



/ o-d/2 r)-2-d/2 

(and a corresponding one when the order parameter is 1 + d/2) to conclude that 



(B.13) 



A s = cP 1 - 1 /« + ... (B.14) 
The medium range solutions in the large z case, Eq. (|3.21|) . are 



l'\,M = p- d/2 { K%W^)' (B.15) 



with [3 = zj{d — 1). The leading order behavior is 



c- 



h ti(ai) 

dh M { ai ) ^cP 1 /^... 

h- M {a x ) =cP- d ^ + ... 

dh M {a x ) = c p-i/e-rf/e + 



(B.16) 



Using these on Am as given by Eq. (|B.6|) . we see that to leading order 
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A M = cP d 't+ 1 ) 

which goes to zero. Therefore we have 



(B.17) 




(B.18) 



Note that, notwithstanding the fractional powers appearing above, is a single valued function, indeed 
near z = it behaves as » z/(d— 1). One also notes that in the large z case A^ is always positive, 
since the Bessel functions I are positive for positive parameter and argument. 

We may perform a similar analysis for the small z approximation, based on Eq. (|3.22|) . 



h% 1 ( P )=p- d ' 2 { ^(7^), (B.19) 



with 7 = 2-v/d/£ + 1. The leading order behavior is 



h M( a i) =c+... 

dh+( ai ) ^cP 1 - 1 ^ + ... 

h M ( ai ) -cP- d /« + ... 

dh- M { ai ) = c p-i/c-rf/e + 



(B.20) 



Using these on Am (cf. Eq. (|B.6fl ^. we see that once again Am behaves at leading order as given in Eq. (|B.17(1 
meaning that it goes to zero as P goes to zero. Therefore we have 



A, 



dh + M {l) 
h + M {l) 



J m (2y/djt + l) 



(B.21) 



In the particular case of £ = 2 the medium range solution can be explicitly calculated for any z, and we 
have 



htt,xiP) = { Y%{^ P )i (B.22) 

where now ft = 2(d + 2) — z/(d — 1). The approximations in Eq. (|B.16(1 or l|B.20|) (for £ = 2 those two 
coincide) are valid unfiromly as £ goes to 2, so when (3 is of order unity, the leading order behaviour 
Eq. I|B.17|) is valid. Note that for z — z^ we indeed have ft of order unity. Now one deduces that A^ = 
— yfftJd/2+i{Vft~)/Jd/2{Vft), and for z = z 2 one finds \fft = (d 2 — d + 4)/4/(d — 1). One can then verify 
numerically that for £ = 2 and z — zi and relevant values of d (between 3 and 8 inclusive) we have A^ > 0. 
By continuity, positivity carries over to values of £ close to 2 and the corresponding dynamo growth rate 
z. This permits us to use near £ = 2 the upper bound w' on the largest solution of Eq. (|4.3|) . and obtain 
the asymptotic behaviour of Sect. 14.21 

B.2 P>1 

Now we have 0,2 = ((d — l)P/2) . The short range solution is in this case 
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with 



hs(p) = P~ d/2 h/2 (VPa'p) 



(B.23) 



a' = -^-=^fz + 2{2- d- d 2 ). (B.24) 
v2 

Similarly to the ?<C 1 case, 

A s = VPa' ^= =cVP+... (B.25) 

The medium range solutions are now power laws, 

h M = p-^ 2 - 2 /^- 1 ^, (B.26) 

where 



Since dhj I {a2) oc V Phf 4(0,2) 



that is, 



y/d(rf3 - io rf 2 + 9 d + 16) + 4(rf -1)1 

5 = 2jd—T) • (R27) 



dh±{a 2 ) - A s h±(a 2 ) = cVPh±{a 2 ) + ..., (B.28) 



A M =c^M + ... = 4 + .... (B.29) 
This goes to zero as P — * 00, and we have 

d-1 2 

where £ was defined in Eq. (|3.16|l . In fact we wouldn't have needed to worry if the limit of Am was infinite 
or zero. The difference would only be a different sign of C, which doesn't affect anything since it is the 
presence of the branch cut alone which determines the positive part of the spectrum. 

C Some Sturm-Liouville theory 

Consider the following general second order linear eigenvalue problem, where a, 6, c are positive functions 
and zeK: 

a[p)h"{p) + b{p)h\p) + c{p)h{p) = zh(p) . (C.l) 
Introduce g = h'/h, then g verifies the first order non-linear (Riccati) differential equation 

/ z-c-bg- ag 2 

g = . (C.2) 

a 

Note that a zero of h corresponds to a pole of g, and the pole is always such that as p increases g goes to 
— 00 and comes back at +00 (since if h is positive before crossing zero then its derivative must be negative 
and vice versa). 
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C.l Montonicity of solutions in z 

Consider for Eq. I|C.1|) the initial condition h'(0) = and h(0) > which in particular implies g(0) = 0. 
Now consider Eqs. (|C.1|1 and (|C.2|) for two different values of z, say Z\ and Z2, and denote the corresponding 
solutions by hi, gi and h%, g 2 respectively. We show that if Z\ > Z2, then gi(p) > g 2 (p) for P less than the 
first zero of h 2 . 

This can be seen as follows. First, the assertion is true near p = since g[(0) — (zi — c(0))/a(0) > 
(z 2 — c(0))/a(0) = g' 2 {0) while gi{0) — 52(0) = 0. Now suppose that at some point the ordering of g\ and 
g 2 changes, this means that the two have to cross, i.e. for some p we have gi(p) = g 2 {p) — G. However at 
this point g[(p) = (z 1 - c(p) - b(p)G - a(p)G 2 ) / 'a(p) > (z 2 - c(p) - b(p)G - a(p)G 2 )/ a(p) = g 2 {p), meaning 
that <?i cannot cross g 2 downwards, which is a contradiction. 

Application 1 

From the above it also follows that the first zero of hi is larger than the first zero of h 2 . Indeed h 2 has no 
zero before its first zero, so g 2 doesn't go to —00 before that point, implying that g\ neither since g\ > g 2 , 
thus hi has no zero either before the first zero of h 2 . 

A particularly useful application of this is to use the position of the first zero of the solution with z = 
as a lower bound on the first zero of any solution for z > 0. 

Application 2 

We may apply the above to the case when Eq. IjC.lf) is taken to be Eq. I|2.15b(l . For P — > the intial 
condition at p — becomes h' M = 0. The case z = can be explicitly solved and we get h]\i(p) = 
p- d / 2 J d ^(2y/d/£ + 1 p?/ 2 ). What needs to be seen is that Hm does not have zeros between and 1, 
equivalent to J d ^ not having zeros between and 2^J d/£ + 1, which follows from the fact that j v> \ > 
2\/v + 1 for v > (where j„ i is the first positive zero of the Bessel function of index u), which may be 
easily verified by a plot or by more serious analysis. We remark that h' M does have zeros, as an effect of the 
term £(d — l)(d + ClP^ ^^M in Eq. I|2.15bj) . and thus the fact that is indeed monotonously increasing 
is surprisingly non trivial because of this term. 

All this allows us to conclude that A^ = /i^ v/ (1)//im(1) grows with z (equivalently with w). 

Application 3 

Along the same lines one can prove for our case the standard lore of Sturm-Liouville theory that if the zero 
mode (2 = solution) has no zeros, then there is no eigenfunction with z > 0. 

The idea is that while the zero mode decays near infinity as a power law, any eigenfunction h\ for 
z > has to decay exponentially, so it will be below the zero mode. On the other hand, from Eq. I|C.1|) on 
deduces that h"(0) grows with z, so that hi has to be larger than the zero mode near p = 0. This would 
imply that the two have to cross in the sense that h 2 comes from above and goes below the zero mode, but 
at the crossing point gi would be less than that of the zero mode, which contradicts the above said. 

C.2 Consequences for modified Bessel function 

We wish to prove here that, for pure imaginary A, the slope of pK' x {p) / K\(p) is bounded from above by 
-1 for all p > 0. 

Using notation from the previous subsections, introduce f(p) = pg{p). Then Eq. l|C.2p translates to 
/' = \{z — c)p + (a/ p — b)f — (a/ p)f 2 ]/a. Applied to the particular case of the modified Bessel equation 
with parameter A 

p 2 h" + ph' - p 2 h = \ 2 h, 
i.e. when a(p) = p 2 , b(p) = p, c(p) = —p 2 and z = A 2 , we obtain 

r = (p 2 + x 2 -f 2 )/p. (c.3) 
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Solving Eq. JH3 for /' = -1 gives f 2 = s(p) 2 where we define s(p) = -[{p + 1/2) 2 + A 2 - 1/4] 1 / 2 . 
Moreover when f 2 > s 2 then /' < — 1. 

We now take / = pK' x (p)/ K\(p) in the case when A is pure imaginary. Then, for large p, asymptotically 
f(p) — s(p) ~ —(1 — 4A 2 )/(16p 2 ) < 0, the last inequality being guaranteed by the fact that we consider the 
case when A is pure imaginary and hence A 2 < 0. This means that for large p asymptotically f < s. 

Using the fact that / is continuous, if / were to become larger than s for some finite p, necessarily it 
would pass through f = s, but at that point we would have f = —1 > s' (the inequality holding for A pure 
imaginary), which is a contradiction to the fact that for larger p we should have / < s. 

This proves that / < s < when s is real, and thus f 2 > s 2 for all p > 0, whence /' < —1 for p > 0. 

C.3 Real spectrum 

Though we do not consider M to be self-adjoint, its spectrum is always real, for the following reason. 

Since M. is a second order differential operator we may conjugate it by a multiplication operator (by 
a "function" which is known in the thory of diffusion processes as the speed measure) to get a symmetric 
operator M. , and taking into account the boundary conditions we have (we see that for any zgC \ R_ 
the solution of Aih = zh which verifies the boundary conditions is a twice differentiable function with zero 
derivative at p — and exponentially decaying as p — > oo, so h is also in the domain of M.'), we can use 
the same trick as for self-adjoint operators: suppose Mh = zh and write J hMh = z J hh, now take the 
complex conjugate of both sides, and since M is real and symmetric, we have J hMh — z J hh, showing 
that z = z, i.e. that z is real. 



D Exact results for P = 

Here we want to study more rigorously the case of P = 0. In this case we can find exactly the zero mode 
of Eq. H2.7JI and show that for A real (recall its definition from Eq. H3.12(l ^ it has no nodes. On the other 
hand for A pure imaginary it has an infinity of nodes. 

Recalling Eq. (|2.11|) . first we have to solve the zero mode equation [£(d — l)(d + £,)p^~ 2 + (d 2 — 1 + 
2£)pt- 1 d p + (d- l)p^d 2 } + [(d 2 - l)p~ 1 d p + (d- l)d 2 }. At Prandtl number the boundary condition is to 
have finite limit at p = 0. The appropriate solution is 



(p« + i)(«*-3)/(«*-D 2 iPi 



2(d-2)4+d(rf-l) + (2-.;)(d-l)A 2(d-2)g+d(d-l)-(2-g)(d-l)A 
H(d-1) ' 2{(d-l) 



where 2F1 is the hypergeometric function, which we shall simply denote as 2^1(0, b; c; a;), and a, 6, c, x are 
defined accordingly to the above displayed formula. 

Let us start with the case of A real. Without loss of generality, we may suppose A > (or otherwise 
exchange a and 6, since the hypergeometric function is symmetric in those arguments). Notice that 2(d — 
2)£ + d(d - 1) > and [2(d - 2)f + d(d - l)] 2 - [(2 - £)(d - I)A] 2 = 2d 2 {d - l) 2 - 8£ 2 (d - 2) > 2[d 2 (d - 
l) 2 - 16(d - 2)] > 0, implying b > 0. 

Notice also 2(d - l)(d + f) - [2(d - 2)£ + d(d - 1)] = d(d - I) + 2£ > implying c - b > 0. 

Now write the following integral representation of the hypergeometric function: 

whence 2-^1(0, b] c; x) > for any x < 1 and c > b > 0. 

Since we have shown above c > b > and since our x < 0, this proves that the zero mode has no zeros, 
and hence there cannot be a dynamo effect. 

For A pure imaginary it is possible to make a large p development using the so called linear transformation 
formula 2 i ;i i(a,fo;c; -x) = p|^p^°| x~ a 2 J 1 i(a, 1-c+a; 1-b+a; -I / x) + ^^^- ) x- h 2 F 1 (b , l-c+6; l-a + 
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b;—l/x). Since a and b are complex conjugates in the case of pure imaginary A, the large x asymptotics can 
be written as 2 Fi(a,b; c; —x) ~ ^( r\b)r\c-'a) x ~ a )> which has an infinity of zeros since a has an imaginary 
part. 
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